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When thermal QCD crosses the critical temperature from below, its pressure den- 
sity rises drastically, consistent with the picture of deconfinement and the release 
of partons as light degrees of freedom. On the other hand, the concept of partons 
^— ^ , is a perturbative one, whereas interactions with the infrared modes in the plasma 

■-.^ ■ always introduce non-perturbative contributions. Here I show how partonic cor- 

' relators can be defined in a gauge invariant and non-perturbative manner which 

Cl^' applies to all energy scales. In particular, I compute the magnetic mass for hot 

SU(2) gauge theory and find m,A = 0.36(2)g 2 T, whose inverse for large T is the 
largest correlation length in the system. 



1. Introduction 



After many years of theoretical studies and well into the period of heavy ion 
collision experiments trying to establish the quark gluon plasma of QCD, we 
are still far from understanding which objects constitute the fundamental 
degrees of freedom in that phase, and what their properties are. Lattice 
simulations of the equation of state have provided firm evidence that the 
effective number of light degrees of freedom is growing rapidly across T c ~ 
170 MeV, as shown in Figure 1. This conclusion can be drawn because the 
pressure rises with temperature for a theory with a given particle content, 
as well as with the number of light fermion flavors for all temperatures. In 
the naive deconfinement picture this is explained by the hadronic degrees 
of freedom dissolving into partons, for which one expects larger correlation 
lengths. 

While this picture explains the qualitative features, the flattening of the 
pressure short of the ideal gas value indicates that up to T ~ 4T C inter- 
actions still play an important role. This in accord with the fact that the 
running gauge coupling at these temperatures is still large, and confirmed 



1 



February 1, 2008 20:16 WSPC/Trim Size: 9in x 6in for Proceedings 



scwm 



2 



5 







4 



3 



2 




100 



200 



300 



400 



500 



600 



Figure 1. Temperature and flavor dependence of the QCD pressure density . 

by detailed studies of screening masses in the plasma, which behave non- 
perturbatively: contributions from the soft magnetic modes ~ g 2 T over- 
power those of the electric modes ~ gT at the temperatures in question 2 . 
Eventually the ideal gas pressure is obtained at asymptotically high tem- 
peratures , and the screening masses are dominated by the perturbative 
contributions. But even in this regime the soft modes ~ g 2 T are, through 
dimensional reduction 4 , described by a 3d confining theory whose pertur- 
bation theory is infrared divergent (Linde problem) 5 . So we have evidence 
for the effective degrees of freedom getting more and lighter corresponding 
to some kind of constituent, but because of the interplay between soft and 
hard modes a purely perturbative parton picture is not appropriate. 

To interpolate smoothly between hadronic and partonic regimes thus re- 
quires a non-perturbative study of the dynamics of color degrees of freedom 
which is valid for all scales. Color dynamics is encoded in Green functions 
of quarks and gluons, which in general are not gauge invariant. In perturba- 
tion theory one fixes a gauge and studies e.g. the field propagators directly. 
While these are not physical observables, they nevertheless carry physical 
information about the parton dynamics in their singularity structure. For 
example, the pole mass defined from the quark propagator is gauge inde- 
pendent and infrared finite to every finite order in perturbation theory 6 . 
A similar result holds for the gluon propagator, provided an appropriate 
resummation of infrared sensitive diagrams has been performed 7 . This has 
been used to define the Debye mass mjj ~ gT in analogy to QED from the 
electric gauge field propagator 8 . Similarly, gauge invariant resummation 
schemes have been designed to self-consistently compute the pole of the 
gluon propagator in three dimensions 9 ~ 12 , corresponding to a "magnetic 
mass" regulating the non-abelian thermal infrared problem. However, it 



February 1, 2008 20:16 WSPC/Trim Size: 9in x 6in for Proceedings 



scwm 



3 

has so far not been clear whether these poles exist non-perturbativcly. 

On the other hand, on the lattice the study of partonic Green functions 
is hampered by several problems. It is difficult to fix a gauge uniquely 
and avoid the problem of Gribov copies 13 . Moreover, most complete gauge 
fixings (e.g. the Landau gauge) violate the positivity of the transfer matrix, 
thus obstructing a quantum mechanical interpretation of the results. For 
these reasons it has been argued to focus on correlators of local singlet 
operators only. A non-perturbative definition of the Debye mass in terms 
of singlet operators has been given 14 , identifying it as the lowest screening 
mass in a singlet channel odd under Euclidean time reflection. The infrared 
cut-off on the magnetic scale g 2 T would in this picture be given by a 3d 
glucball mass. However, the corresponding correlation lengths arc hadronic 
and not directly related to screening phenomena like e.g. J/i/'-suppression, 
which are caused by charged intermediate states 15 . 

In this contribution I want to demonstrate that non-perturbative and 
gauge invariant information is contained in an appropriately defined gluon 
two-point function, thus permitting to arrive at a field theoretical defini- 
tion of an associated correlation length 16 < 17 . Sections 2-4 show how to 
construct such an object which can be proved to decay with eigenvalues of 
the Hamiltonian. Section 5 relates the lowest eigenvalue to a level splitting 
between static mesons, which can be used to compute it without recourse 
to any gauge fixing at all. Such a computation is presented for SU(2) gauge 
theory in 3d. In Section 6 the result is compared to those from analytic re- 
summation schemes, and argued to constitute the largest correlation length 
in the thermal system. Section 7 presents the conclusions. 

2. A non-local gluon operator 

A gauge invariant lattice gluon correlator can be defined when a complex 
N-plet transforming in the fundamental representation is available. One 
possibility is to take the eigenfunctions of the spatial covariant Laplacian, 
which is a hermitian operator with a positive spectrum, 

-{Dl[U\) a p$\x) = \ n f£\x), A">0. (1) 

They provide a unique mapping U — > f[U] except when eigenvalues are 
degenerate or |/| = 0. In simulations the probability of generating such 
configurations is essentially zero. These properties have been used previ- 
ously for gauge fixing without Gribov copies and to construct blockspins 
for the derivation of effective theories 18 . The lowest eigenvectors are used 
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(f^\x),f^ 2 \x)), which transforms 



to construct the matrix Q(x) = ttz-t 
as f2 9 (x) = g(x)Q(x)h' f (t). Composite link and gluon fields are defined by 

V„{x) = n\x)U^(x)n(x + fl), (2) 

^T r (V,(x)-V}(x)) 



i 

25 



V„(x)-V}(x)--Tr (V,(x) 



(3) 



both transforming as Of(x) = h(t)Oi(x)h^ (t), whereas Vq(x) = 
h{t)Vo{x)h) (t +1). The gauge field at a given time now has only a global 
gauge freedom left. To cancel this in the correlator, the zero momentum 
projected time links V a (t) = J2 X ^o( x > t)/\ ^ x ^o( x > *)l are multiplied to 
"strings" Vo(t\,t2) connecting two timeslices. These ingredients can be 
combined to the gauge invariant operator 

0[U] = Tr [A 4 (x,0)V/ (0,t)A 4 (x,i)t/ t (0,t)] , (4) 

which in the particular gauge Vo(t) = 1 reduces to the gluon propagator. 
Hence we have a gauge invariant non-local observable, which is identical to 
the gluon propagator in a particular gauge. 

Using the transfer matrix formalism 19 , this observable can be converted 
to a trace over quantum mechanical states, falling off exponentially as 

-E )t 



<O[[/])=]T|(0|i Q/3 (x)|n 



L0\|2 e -(B„ 



(5) 



provided that there is a finite energy gap (Ex — E ). The eigenvalues E n 
and eigenvectors \n L0 ) are those of the Hamiltonian in Laplacian temporal 
gauge, Vo(x) = 1. It has been proved 16 that this Hamiltonian has the same 
spectrum as the Kogut-Susskind Hamiltonian Hq 20 , which is obtained by 
quantizing the theory in temporal gauge. Note that the zero momentum 
projection of the Vq switches off the divergent self-energy of the sources, so 
the energies are finite. 



3. Discussion and interpretation 

One may now ask to what extent these results depend on the particular 
choice of fi[I7], which is not unique. It is crucial that Q depends only 
on spatial links to preserve the transfer matrix. Clearly, any £ SU(N) 
local in time and transforming in the same way permits construction of 
the gauge invariant observable Eq. (4). From the spectral representation 
it follows that all such observables fall off with the same spectrum, while 
f2 only enters the matrix elements representing the overlap of the operator 
with the eigenstates. 
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The construction of the composite link variable Eq. (2) may also be 
viewed as fixing Laplacian gauge on each timeslice. Our result then im- 
plies that all gauges satisfying the mentioned constraints (e.g. the standard 
lattice Coulomb gauge) will produce the same exponential decay. Fourier 
transforming Eq. (5) to momentum space, we obtain the non-perturbative 
analogue to the situation in perturbation theory: the energies appear as 
gauge invariant poles, while the matrix elements correspond to gauge de- 
pendent residues. 

Note that this does not imply an asymptotically free colored state. 
The temporal links in the operator indicate the presence of static charges, 
as e.g. in a Wilson loop. In a Hamiltonian formulation 20 ' 19 the Kogut- 
Susskind Hamiltonian H n acts on a Hilbert space H of all complex, square 
integrable wave functions. Gauge invariant wave functions, ip{U 9 ] = tp[U], 
form a subspace H C H 0} on which the gauge invariant projected Hamil- 
tonian H acts. Generally, the spectrum of H consists of the physical par- 
ticle states of the theory, which couple to local gauge invariant operators. 
In addition to these states, Hq contains also the spectrum of gauge field 
excitations in the presence of static sources, such as the static potential, 
gluelumps etc. These are gauge invariant energy eigenvalues to non-trivially 
transforming pure gauge wave functions. Full gauge invariance is restored 
once the source fields are made explicit. Our correlator thus probes a gauge 
field excitation with the quantum numbers of the gluon in the presence of 
sources, which do not contribute to the field energy. 



4. A numerical check 

It is expedient to first test this new operator in a 3d SU(2) Higgs model, 
which in its broken phase has physical states with the quantum numbers 
of the gluon, the W-bosons, with detailed results available 21 . Figure 2 
compares the W-mass measured in these works by the standard operator 
V[U, <f>] — Im (x)U(x)4>(x + fi)), with results obtained from the com- 
posite links T^[C/]. Full agreement is observed for different lattice spacings, 
which also confirms that the energy levels Eq. (5) do have a continuum 
limit. On the other hand, when simulated in the pure gauge theory, they 
grow linearly with the box size 17 . This effect stems from the no n- locality 
of the functional £l[U], which depends on all link variables in a t-slice. On a 
periodic torus and in a confining regime, it will thus project predominantly 
on torelonic states and be blind to local quantities. 
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Figure 2. The W-boson mass in a 3d Higgs phase, computed from the standard 
operator V[<j>, U] 21 and the non-local V[J7], Eq. (2) 17 . The gauge coupling is 3d. 



5. The gluon propagator and static mesons 

It is then necessary to construct an alternative operator coupling to the 
states we are interested in. This is achieved by introducing explicit scalar 
fields for the fundamental static sources. In continuum notation, adding 

S4u,<f>] = Yl {-\ D M*)\ 2 + ™%\<K*)\ 2 } (6) 

to the pure gauge theory results in QCD with scalar quarks. In the 
limit m — > oo the scalars become static sources propagating in time 
only, their propagator being known exactly to consist of temporal Wil- 
son lines f7o(*i)*2)- Scalar and vector mesons are described by S(x) — 
(f>^(x)<p(x),V(x) = lm((f)^(x)Di(x)<p(x)), respectively. In the static limit 
the scalar fields do not contribute to angular momentum, nor can they be 
excited into higher quantum states since they are quenched. Consequently 
the mass difference vtia = lim m0 ^ oo (My — Ms) is a pure gauge quan- 
tity, characterizing an excitation with the quantum numbers of the gluon. 
Moreover, integrating the scalars out analytically, one obtains for the ratio 
of correlators 

(V( X ,0)V(x,t)) c I DATl (A(x,0)^ (0,t)A(x,t)t4(0,t))e-^« 
(S(x,0)S(x,i)) c JDATi (Uo(P,t)U%(0, t)) e~ s ™ 

(7) 

In temporal gauge this reduces to the gluon propagator again. In other 
words, the mass difference My — Ms in the static limit should be equivalent 
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Figure 3. The lowest states in 3d scalar QCD 1 . Mq denotes scalar glueballs, Mg y 
scalar and vector mesons. The gauge coupling is 3d. 



to the pole mass of the gluon. (Similarly, a correlation length for thermal 
electric gluons was recently extracted from a ratio of singlet correlators 22 .) 

Numerical results for 2+1 dimensional SU(2) are shown in Figure 3. 
With increasing scalar mass, the measured glueball states Mq attain their 
pure gauge values indicated by the dashed lines. The scalar bound states 
move out of the spectrum, with My — Ms approximately constant. At the 
largest scalar mass one finds itia — 0.37(6)g 2 , or MqItua ~ 4.2. 

Rather than approaching the static limit numerically, one may also 
compute it directly, beginning from the right hand side of Eq. (7). The 
numerator is a continuum version of the lattice operator Eq. (4) , the dif- 
ference being that the time-like strings are not zero momentum projected. 
In discretizing the covariant derivative, the non-local functionals should 
be avoided. This is achieved by observing that the exponential decay of 
a correlator is entirely determined by the Hamiltonian and the quantum 
numbers of the operators used. Instead of discretizing the covariant deriva- 
tive, one can then equally well employ a local higher dimension operator 
sharing the same quantum numbers, such as (DjFij(x)) a . The discretized 
version of the latter is simply the adjoint part of a linear combination of pla- 
quettes, and hence a local operator. The denominator, on the other hand, 
consists of a Wilson line running back and forth, i.e. an adjoint Wilson line, 
which has the same exponential decay as a field strength correlator 23 . The 
discretized version of Eq. (7) is then simply given by a ratio of gluelump 
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correlators with the appropriate quantum numbers, 

(lr (A^(x, 0)(7 (0, t)A^(x, t)t4(0, t)) ) 

(Tr (fi J -(x,0)I7 (0,t)fi J -(x, t)t4(0, t))) 

and its asymptotic exponential decay is the mass difference of the corre- 
sponding gluclump masses, in which the divergent self-energies of the Wil- 
son lines cancel. (Note that, had we coupled the gluon to one adjoint source 
field instead of two fundamental, we would have obtained the same static 
limit.) The calculation performed in 17 and extrapolated to the continuum 
gives as the final result 

m A = 0.360(19)^. (9) 
Note the agreement with the result obtained by using heavy scalar fields. 

6. The magnetic mass 

In the framework of the 4d thermal field theory, the 3d gluon mass appears 
at asymptotic temperatures as the magnetic mass for hot SU(2) gauge 
theory, which then is tua = 0.36(2)g 2 T. It is now interesting to compare 
with gauge invariant resummation schemes of perturbation theory and a 
Hamiltonian strong coupling analysis, all done in the 3d gauge theory, which 
have been used in the past to compute the pole of the propagator. As the 
Table shows, the leading order results of all these calculations get within 
20% of the right answer, and the two-loop calculation 24 in one of the 
schemes 9 even suggests reasonable convergence. 





Ref. 


m A/gl 


1-loop gap eq. 


10 


0.38 




9,11 


0.28 




12 


0.25 


2-loop gap eq. for 9 


24 


0.34 


Hamiltonian strong coupl. 


25 


0.32 



For non-pcrturbative evidence for the role of this quantity in non-abelian 
plasmas, recall Guy Moore's contribution to SEWM 20 26 , in which he 
showed the finite size scaling behavior of the sphalcron rate in the hot SU(2) 
pure gauge theory. This is obtained by simulations of the effective theory 
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Figure 4. Volume dependence of the sphaleron rate, which picks up at 2l ma g (the 2 is 
for the periodic boundary conditions) 26 . 



for the soft modes ~ g 2 T and shown in Figure 4. A quantity is typically 
afflicted by finite size effets as long as the correlation length corresponding 
to the relevant modes does not fit into the simulated box, L < £, and begins 
to approach its infinite volume limit once L > £. The sphaleron rate dis- 
plays a rather clear signal for this, and the correlation length one estimates 
from the plot is fully compatible with rn^ 1 , while being completely at odds 
with Mg 1 . 



7. Conclusions 

By constructing a non-local lattice operator whose correlation function 
is amenable to the transfer matrix formalism, 1 have shown that a non- 
perturbative, gauge invariant mass scale is associated with the gluon field, 
representing the smallest eigenvalue of the Kogut-Susskind Hamiltonian in 
the presence of external charges. In momentum space it corresponds to 
a pole in the propagator in all unique gauges that are local in time. The 
eigenvalue can be related to some particular level splittings of static mesons. 
I have calculated this quantity for the 2+1 dimensional SU(2) pure gauge 
theory, and found it to be roughly a quarter of the glueball mass. This 
implies a magnetic mass for the hot SU(2) gauge theory, whose inverse 
constitutes the largest correlation length of the system at asymptotic tem- 
peratures. Of course, close to T c the increase of g 2 (T) may lead to level 
crossings with electric modes and change this picture, and so will the ad- 
dition of light fcrmion flavors. But non-perturbatively dressed partons are 
theoretically accessible and, as the sphaleron example shows, might be the 
relevant degrees of freedom for certain aspects of plasma dynamics. 
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